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Abstract 

The rapid development of micro- and nanooscillators in the past decade has led to the emergence of novel 
sensors that are opening new frontiers in both applied and fundamental science. The potential of these novel 
devices is, however, strongly limited by their increased sensitivity to external perturbations. We report a non- 
invasive optomechanical nano-stabilization technique and apply the method to stabilize a linear nanomechanical 
beam at its ultimate thermodynamic limit at room temperature. The reported ability to stabilize a mechanical 
oscillator to the thermodynamic limit can be extended to a variety of systems and increases the sensitivity range 
i—l of nanosensors in both fundamental and applied studies. 

43 Thermal noise is well known to be a ubiquitous and fundamental limit for the sensitivity of micro- and nanome- 

chanical oscillators [U El E] . This noise arises from the unavoidable coupling of the mechanical device to a thermal 
bath [3] , which is responsible for an incoherent motion with random amplitude and phase 5 . A number of precision 
£h measurements based on using nanomechanical resonators, such as gradient force [5J[7], mass [SUS] and charge [TU] 
, t detection, or time and frequency control |llj . consist in detecting linear changes induced by the measured system 
on the phase of the coherently driven nanomechanical sensor p] [12]. These applications are thus theoretically 
limited by the random phase fluctuations related to thermal motion. Decreasing the effect of thermal noise can 
be accomplished by increasing the relative contribution of the coherent drive pQ. However, the latter has to be 
gH kept below a certain threshold, above which the mechanical oscillator becomes nonlinear [3J [UJ : At this threshold, 
q thermal noise constitutes an irreducible sensing limit for linear nanomechanical measurements. In practice, even 
O when driven close to their nonlinear onset, nanomechanical resonators presently operate far away from the thermal 
limit [H HH [TU [16] . This is due to a plethora of external perturbation sources acting on the oscillator E| , such 
as temperature fluctuations [T7], adsorption-desorption noise [T5] or molecular diffusion along the oscillator [TS] , 
J> While several proposals have been made [2"0l |2"T] and implemented Ell E2] to decrease excess phase noise, they are 
based on oscillators being driven into the nonlinear regime and require specific operating conditions, thus making 
the schemes difficult to adapt for a general class of resonators. In contrast, here we introduce an approach that 
generally applies to linear nanomechanical resonators. This approach employs an auxiliary mechanical mode as a 
• frequency noise detector, whose output is used to stabilize the frequency of the fundamental out-of-plane mode of a 
nanomechanical beam. Its experimental implementation enables quasi-optimal linear operation at room tempera- 
ture, with an almost perfect external phase noise cancellation, as demonstrated for the nanomechanical beam driven 
i-H just below its nonlinear onset. Hence, our system operates close to its ultimate thermodynamic linear sensing limit. 
IL^ Our scheme can be readily incorporated into a variety of other systems used for applications in [2l [9] , frequency 
control [TT] or fundamental studies EH ES] . 

Our system consists of an integrated hybrid optomechanical transducer described elsewhere |26j . Briefly, it 
consists of a 90 /im x 700 nm x 100 nm low dissipation nanomechanical beam made of high-stress silicon nitride 
(SiaN4) placed in the near-field of a high-Q whispering gallery mode (WGM) confined in a disk-shaped microcavity 
(cf. Fig. 1(a)). Laser light is coupled in and out of the oscillator with a tapered fiber. Mechanical motion is read 
out using the optomechanical interaction by placing the laser frequency at mid-transmission of the optical mode, 
such that the phase fluctuations of light induced by the mechanical motion are directly transformed into amplitude 
fluctuations readily detected with a photodiode EZ] . 
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Our stabilization scheme is based on the intuition that most of the observed external frequency noise in our 
system arises from perturbations affecting the entire beam geometry (e.g. temperature fluctuations causing causing 
changes in length, strain among others). Under such conditions, the induced frequency fluctuations of the various 
cigenmodes are expected to be highly correlated. One of these modes can thus be used as a "noise detector" 
whose output serves as an error signal in a feedback mechanism applied to the beam, enabling a cancellation 
of the geometry variations and thereby a stabilization of all eigenmodes (cf. Fig. 2(a)). In the following our 
analysis will focus more specifically on two modes, which are the lowest-order out-of-plane mode, which we seek 
to stabilize, and the lowest-order in-plane mode, which is used as the noise detector. These two modes (referred 
to as mode 1 and 2, respectively in the following) have resonance frequencies vm,i = ^m.i/27t = 2.84MHz and 
^m,2 = ^m,2/27t = 3. 1 05 MHz as well as linewidths of rM.i/27r = 7.5 Hz and Tm^/Zk = 48 Hz. Note that using the 
in-plane mode as the noise detector has several advantages: First, the accurate position resolution of the beam's 
intrinsic motion enables to achieve a large bandwidth (i/rj ~ 20 x Fm,2/27t ~ f kHz, cf. Fig. 1(d)) over which 
frequency fluctuations can be detected with thermal limited imprecision. Second, this geometry enables a selective 
exposure of the perpendicular out-of-plane mode to external signals, thus preserving its sensing capabilities when 
feedback is turned on. Experimentally, the frequency fluctuations of each mode are detected via its out-of-phase 
response to an external radiation pressure coherent drive with a frequency of Qp,, (i € {1, 2}), similar to a phase- 
locked-loop detection. We choose the driving strengths such that both modes have an amplitude of just below the 
respective nonlinear thresholds. Both measured out-of-phase quadratures for mode 1 and 2 are shown in Figure 
2(d) with the demodulated response plotted versus the detuning LUi/2ir = (flpj — ilM,i)/2 7r between the driving 
and the resonance frequency. By putting the external drive close to the resonance frequency (Qp^ ~ ^m,i)j we 
can monitor the resonance frequency as small changes on the order of Sfli < Fjy^; that are proportional to the 
out-of-phase quadrature signals. Our feedback mechanism relies on DC-photothermal actuation applied to the 
beam by sending a third amplitude modulated laser tone (control beam) with input power in the mW range (see 
Fig. 1(e)). The frequency shift follows linearly the optical input power as shown in Figure 2(b) for mode 1. 

Figure 2(c) shows the frequency response of modes 1 and 2 to a triangular wave function applied to the control 
beam. One can already note the similarities of the fluctuations superimposed to the response triangles in both 
cases, indicating the common origin of the frequency noise over the various eigenmodes. 

These noise correlations can be quantitatively confirmed by determining the frequency correlations Ci_2[w] = 

1/2 

(Svi [w]<5^2 M)/ (Si/Jwj^H) (with fo/jfcd] denoting the Fourier transform of the frequency fluctuations of mode 
i, S Vi its spectrum and (...) statistical average), which are represented in Figure 3(c) (inset). It should be emphasized 
that in absence of external frequency noise, the frequency thermal imprecisions over modes 1 and 2 are expected to 
be uncorrelated: the high level of correlations is therefore a signature of the presence of external frequency noise. 
Experimentally, this excess of frequency noise is determined by comparing the phase noise spectrum of mode 
i to its thermal equivalent imprecision, which in the strong drive limit is given by [2J [5] 



oth r l 1 Tm^bT . . 



with (X 2 SC ■) as the variance of the driven amplitude, m e s^ as the effective mass, &b as Boltzmann's constant, 
T as the temperature and lu as the Fourier frequency. Equation [I] simply reflects the fact that phase noise can be 
interpreted as the " angle" under which the thermal distribution appears from the center of the phase space (see Fig. 
[3](b)) . The time evolution of the frequency drift Svms, as determined from the mechanical out-of-phase response, 
is shown in Figure 3(a) for mode 1 and mode 2. It is clearly evident that the frequency noise is strongly decreased 
for the stabilized case. These time traces are used to obtain the phase noise spectrum via a Fourier transform. 
The results are shown in Figure 3(c), both in the non-stabilized (red) and stabilized (blue) cases. Dark colored 
curves correspond to mode 1 and light colored to mode 2. The results are calibrated using the inferred driven 
root-mean-square amplitudes (^osci) 5 = 1-7 nm and (^osc 2) 5 = 4.1 nm, just below their respective nonlinear 
thresholds. When the system is not stabilized, the phase noise of both modes is clearly above the thermal limit. 
With the stabilization scheme turned on, the phase noise of mode 1 is reduced to the thermal limit in the frequency 
range of ^10-50 Hz. The spectrum of mode 2 implies that its phase noise was reduced to even below the thermal 
limit, which is clearly unphysical. This artifact arises from squashing, which is generally observed in an in- loop 
measurement where the feedback gain is so large that the reinjected background - in our case the thermally induced 
frequency fluctuations - is the dominating part of the error signal. As the spectrum of mode 1 is obtained in an 
out-of-loop measurement, reinjection of background would manifest itself correctly as an increase in the phase noise 
spectrum. Since the thermal frequency noise background — u} 2 S^ i is smaller for mode 2 (28 mHz/vTiz) in 

comparison to mode 1 (41 mHz/\/Hz), the reinjection of noise is limited for mode 1 - the mode of interest - to a 
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level of 1.5 dB with respect to the thermal limit. 

We also observe the effect of stabilization on the correlations between the frequency noise of both modes. As 
noted above, noise suppression is expected to correspond to a significant decrease of the frequency noise corre- 
lations, as the noises resulting from fundamental thermodynamic fluctuations are uncorrelated. We confirm this 
experimentally as shown in the inset in Figure 3(c). 

To access the long-term stability for mode 1, we determine its Allan deviation. In order to avoid measurement 
artifacts, the mechanical motion is no longer driven using an external frequency reference, but is maintained in self- 
oscillation using a feedback circuit [28l [26] (see Fig. 1(e)). The readout signal is fed into a frequency counter (Agi- 

r ,1/2 
lent 53230), which determines the Allan deviation as (<5^m,i/ '^m,i)t = 2(n~i) S/c=i ((^M,fe+i — ^M,k) I v m) > 
with &M.k as the averaged frequency in the kth discrete time interval over the gate time r. 

The obtained results are shown in Figure 3(d), where the red and blue dots represent the values of the fractional 
Allan deviation of mode 1 as a function of gate time in the non-stabilized and stabilized cases, respectively. The slow 
fluctuations are significantly reduced in presence of stabilization, with the noise suppression even exceeding 20 dB 
at the second scale. The achieved long-term stability is further compared to the theoretical limits of our system. 
For a self-sustained mechanical resonator, the thermo-mechanically limited Allan variance <Ji{t) = (Sv^ 1 /^'m,i)t 
can be shown to be given by (see Supplementary section S2.2): 

/ \ ( (^th,l) 1 \ 

The data of Figure 3(d) were acquired in presence of a self-sustained displacement (X 2 SC 1 ) 1 / 2 ~ 100(X 2 h 1 ) 1 ^ 2 . 

This analysis shows that our stabilization scheme enables reducing external frequency noise close to the thermal 
limit, but not closer than a factor of 2 (that is 6 dB in power), which disagrees with the spectral data of Figure 3(c). 
This can be explained by a reduced sensitivity of noise detection: Mode 2 being externally driven, the expression 
of its thermal limited Allan deviation differs from mode 1 and is described by (see Supplementary section 2.2): 



ct 2 (t) 



(*th, 2 ) 3 + e - r M,2T _ 4 g- 



(A Q 2 SC , 2 ) (O m , 2 t) 



(3) 



With (^osc^) 1 ^ 2 — 10 2 (^t 2 h 2) 1 ^ 2 ' the noise detection imprecision (5v^_ 2 )t~o.i a — 10 8 ^m,2 is comparable to 
(^Mi)^o.is — 10~ 8 ^m.i (Eq. [2J), and therefore has to be taken into account when determining the stabilization 
limit, which is given by: 



- («,i> 2 + (-5<2) 2 ) 1/2 - (4) 

The thermal limit for mode 2 (Sv^ 2 )t as well as the overall stabilization limit given by Eq. Hare represented 
in Figure 3(d), showing very good agreement with the experimental data up to gate times in the TOOms range. On 
longer timescales, the frequency fluctuations are stabilized to a flat level in the lOppb range, which we attribute 
to insufficient gain. 

Finally, we analyze the frequency stability and its improvement with the stabilization scheme by using a third 
independent method based on a recent proposal [35] . Briefly, this method is based on the simultaneous measurement 
of both in-phase and out-of-phase quadratures Xi(t) and X2(t) of the coherently driven mechanical motion. These 
quadratures are used to reconstruct the complex motion trajectory u(t) = Xi(t)+iX2(t) whose moments (u n ) can be 
shown to be insensitive to thermal noise, enabling the extraction and quantitative study of external frequency noise. 
Thus, the presence of external frequency fluctuations <5S^m (t) is revealed by an anomaly of the ratio r — (u 2 )/{u} 2 , 
which falls below unity by an amount proportional to the frequency noise power. Indeed, the quadratures of the 
coherently driven motion write as the sum of a steady and a fluctuating part Xi^{t) = (^1,2) +5Xi^(t)- A straight 
forward expansion gives (u 2 ) = (u) 2 + ((SX1) 2 ) — ((SX2) 2 ) — 21(8X18X2) ■ With the thermal contribution to the 
quadrature fluctuations 8Xl\ being uncorrelated and with identical variances, the previous expression becomes 
(u 2 ) = (u) 2 + ((SX^) 2 ) - ((8X^) 2 } - 2i(8X1 xt SX 2 °* t ), with Xf$ denoting the quadratures fluctuations due to the 
presence of external frequency noise. For a resonantly driven system, the in-phase quadrature Xi is maximum, and 
SX1 xt = (dXx/du)) 8n M (t) = 0, (u 2 ) is thus reduced to the quantity ((SX^f) ~ (dX 2 /duj) 2 (SD^) ~ ^fS^(u) 2 - 
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At the mechanical resonance frequency, r gives a direct access to the frequency noise power, with no calibration 
being required. In a more sophisticated approach [29 , it can be shown that for a weak Gaussian frequency noise, 
the ratio r[u>] varies along the driving frequency as 



r[u] = ^Hrk? (5) 

with S UU1 denoting the frequency noise spectral density. Figure 4 shows the experimental determination of r[ui] 
we performed with our system, both in the non-stabilized (Fig. 4(c)) and in the stabilized cases (Fig. 4(c,d)). 
Figures 4(a) and 4(b) show the realizations of the phase-space trajectory used to determined r[ui ~ 0] in the 
non-stabilized and stabilized configurations respectively. 

As expected, without stabilization we obtain a discernible dip around the mechanical resonance frequency as 
shown in Figure 4(c). The dip depth is about 1 — r[ui = 0] ~ 4 x 1CP 2 with the deviations from the expected 
shape resulting from long-term frequency drifts that change 5ui. Following Eq. [5] such a dip is expected to 
correspond to a frequency noise S vv — S LULU /(27r) 2 = (218mHz/v / Hz) 2 , that is about 15 dB above the thermal limit 
S^w] = ctS^.iM ~ (41 mHz/vTIz) 2 , in very good agreement with the results presented in Fig. 3lc). When 
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the stabilization is switched on, the dip is greatly reduced to a level of 1 — r[cu — 0] ~ 6 x 1CP 4 , corresponding to 
a residual noise S vv — (27mHz/\/Hz) 2 that is 1.5 dB above the expected thermal limit. This value confirms the 
origin of this remaining imprecision as resulting from squashing, which is expected to generate the very same noise 
excess as explained above. 

In conclusion, we have proposed and implemented a scheme allowing a non-invasive stabilization of a mechanical 
mode. This technique proves to be remarkably efficient, as shown by the almost perfect noise cancellation for the 
lowest order out-of-plane mode of a high-Q nanomechanical resonator operated just below its nonlinear threshold, 
which places our system close to the optimum noise limit at room temperature. Our method is very general and 
can be readily applied to a number of ultra-sensitive nanomechanical systems: First, our scheme allows using any 
mechanical mode for detection of frequency noise and we successfully demonstrated noise cancellation using the 
third-order out-of-plane mode. To apply an error signal one needs a DC restoring force, which is widely available 
either via photothcrmal forces in optomechanical systems such as in our work or as electrostatic forces in electrical 
systems [501 |3~TI 152"] . Importantly, our scheme fully preserves the detection capabilities of the stabilized system, 
which we have experimentally verified by showing efficient radiation-pressure induced thermomechanical squeezing 
[33] of the fundamental out-of-plane mode [34] that was impossible to observe without stabilization due to long-term 
frequency drifts. 

The fabrication was carried out at the Center of MicroNanotechnology (CMi) at EPFL. Funding for this work 
was provided by the NCCR Quantum Photonics, the SNF, DARPA Orchid and an ERC. 
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Figure 1: The hybrid optonanomechanical system and experimental setup, (a) Optical micrograph 
showing 4 of the 35 hybrid structures present on the chip used in the experiment, (b) Magnified view of (a) 

1/2 

showing the details of the hybrid structure (false colors), (c) and (d) Thermal displacement spectra S xx [^] 01 the 
fundamental out-of-plane (mode 1) and in-plane (mode 2) modes, respectively. The green curves denote the thermal 
motion of the mechanical modes. The insets show finite-element simulations of the respective mechanical modes 
with the red arrows representing the direction of displacement. The gray area in (d) represents the bandwidth 
Avd in which the stabilization scheme enables noise cancellation at the fundamental thermodynamical limit, (e) 
Experimental setup for the stabilization mechanism, with the readout and actuation of the mechanical oscillator 
being obtained via optical means with external cavity diode lasers (ECDL). The readout signal from the photodiode 
(PD) is split into a demodulation circuit for the directly driven noise detector mode 2 and a circuit for mode 1, 
which is either a directly driven (straight line) or self-driven (gray line) circuit. Two amplitude modulators provide 
the means to drive the mechanical modes (AM2) and to control the mechanical frequency by changing the optical 
power coupled into the cavity (AMI). The frequency fluctuations Svyi^i and <5z^m,2 are recorded either with an 
oscilloscope, when the system is operated in the direct driving mode or <5i^m,i is recorded with a frequency counter 
when the mode is self-driven. 
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Figure 2: Detecting and stabilizing frequency noise with an optonanomechanical system, (a) Scheme 
describing the stabilization principle. Phase fluctuations of the in-plane mode (mode 2) are detected by determining 
its out-of-phase response to a coherent drive and serve as an error signal. By applying a feedback control based 
on these fluctuations, the out-of-plane mode (mode 1) is stabilized, (b) Resonance frequency shift Avm.i of mode 
1 in terms of its mechanical linewidth Tm,i versus input optical power, (c) Resonance frequency shifts Az/m,i and 
Avm,2 of modes 1 and 2, respectively, for a saw-tooth modulation of the input optical power P. The correlation of 
the response of both modes is clearly visible, (d) The out-of-phase mechanical response of mode 2 (upper graph) 
and mode 1 (lower graph) to a coherent drive versus the detuning (f2p — ^m)/27t. 
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Frequency (Hz) Gate time (s) 

Figure 3: Frequency stabilization at the ultimate thermodynamic limit, (a) Time evolution of the fre- 
quency drift Sum for mode 1 (upper graph) and mode 2 (lower graph). The non-stabilized case is given by the 
red curve and the stabilized case by the grey one. (b) Thermal equivalent phase noise in phase-space: Xi and X 2 
are the in-phase and out-of-phase components of mechanical motion, respectively, and form the phase-space. In 
this representation, the phase of mechanical motion identifies to the azimuthal coordinate (j). The phase noise S4> 
simply corresponds to the angle under which the thermal distribution appears from the center, (c) Frequency noise 
excess of modes 1 (dark-colored) and 2 (light-colored) over the Fourier frequency components. The red and blue 
curves show the non-stabilized and stabilized case, respectively. The stabilization scheme reduces the frequency 
noise by over 10 dB up to 90 Hz and allows to stabilize mode 1 to the thermal limit for a wide bandwidth. Inset: 
The correlation of frequency response for both modes for the non-stabilized (red) and the stabilized (blue) case, 
(d) The fractional Allan deviation for the non-stabilized (red) and stabilized (blue) case. The dashed line shows 
the thermal limit associated with mode 1 for a self-driven oscillator. The dot-dashed line shows the thermal limit 
for mode 2, which was externally driven. The full line shows the thermal limit for mode 1 in presence of feedback 
(EqRJ by taking into account the reinjected imprecision from mode 2. 
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Figure 4: Frequency noise self-characterization, (a) A realization of the phase-space trajectory acquired in 
the non-stabilized case when driving the nanomechanical beam at the mechanical resonance frequency (oj ~ 0). 
The acquisition duration is on the order of one minute. The red Gaussian disk represents the thermal distribution 
of variance (X 2 h ), responsible for an equivalent phase noise (5</> t h-(b) Same as in (a), but for the stabilized case. 
The remaining phase noise clearly appears as being substantially reduced compared to non-stabilized case, (c) The 
ratio (u 2 ) / (it) 2 for mode 1, which signifies the presence of frequency fluctuations due to non-thermal contributions, 
versus detuning from resonance frequency (red: nonstabilized; blue: stabilized). The red dashed-line corresponds 
to the model given by Eq. [5] for a Gaussian noise with frequency power spectral density 15 dB above the thermal 
limit. Stabilization greatly reduces the dip due to excessive non-thermal noise, (d) Detailed view of (u 2 )/(u) 2 
for the stabilized case. The blue dasehd-line correspond to the model given by Eq. [5] for a Gaussian noise with 
frequency power spectral density 1.5 dB above the thermal limit. 
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Supplementary Material - Stabilization of a strongly driven nanomechanical 

oscillator to the thermal limit 



1 Phase noise spectrum of an oscillator at thermal equilibrium 

In this section, we present the calculation of a thermal equivalent phase noise for a mechanical resonator at thermal 
equilibrium. Both cases of a free running and externally driven oscillator are treated successively. In the following, 
we assume the oscillator to be well described within the weakly damped harmonic approximation, with mechanical 
susceptibility x[f2] given by: 

1 



XM = 



(SI) 



with M, Qm and Tm as the mass, mechanical resonance frequency and damping rate respectively. 
1.1 Free-running thermally driven oscillator 
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Figure SI: Phase-space representation of a mechanical state at thermal equilibrium, (a) Phase-space trajectory of a 
mechanical oscillator left at thermal equilibrium. The trajectory describes in cartesian coordinates (Xi(t), X 2 (t)), 
or equivalently in polar coordinates (e(i), <p(t)). (b) Phase-space trajectory of a mechanical resonator driven at the 
mechanical resonance. The thermal state is "displaced" of a quantity (^osc) 1 ^ 2 a l° n g the abscise. 

We start this calculation with noting that the phase <p(t) of a thermal state x(t) can be expressed as a function 
of its slowly varying in-phase and out-of-phase motion components Xi(t) and X 2 {t) (x(t) = X\{t) cosJ7m£ + 
X 2 (t) sinf2]y[f, ^m/27t denoting the mechanical resonance frequency) as follows: 

= 2arctan(! 2 H). (S2) 

Here the factor of 2 preceding the arctan function is to allow the phase to vary over the interval [0, 2tt\. In order to 
obtain a suitable expression for deriving the spectral properties of ip, we subsequently take the derivative of Eq |S 2[ 
which gives: 

ie(t)£(t) = X 1 (t)X 2 (t)-X 2 (t)X 1 (t), (S3) 

with e(t) = Xf(t) + X 2 (t) being proportional to the energy operator. The phase and the energy being independent, 
the autocorrelation of the quantity op can be factorized, which leads to: 

±{e(t)e(t')){<p(t)m) = (Mt)XiV))<Mt)M*)) + (X 2 (t)Mt'))(X 1 (t)X 1 (t')) 

-(X^X^iX^X^t')) - (X 1 (t')X 1 (t)){X 2 (t')X 2 (t)). (S4) 
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In the above equation, (...) is for statistical average, t and t' for two arbitrary moments in time, and the factorization 
on the right side arises from the well known statistical independence of the motion quadratures X\ and X 2 . For 
evident reasons of symmetry, the two terms on the right side of Eq. S 4 on the first line (resp. on the second line) 
are identical such that it can be written: 

\W)<*))W)W)) = (Mt)x 1 (t'))(x 1 (t)x 1 (t')). (S5) 



We now determine each factor appearing on the ride side of Eq. S5 We start with the simpler one (X 1 (t)X 1 (t')), 
which identifies with the autocorrelation function of the quadrature X\. Assuming the stationarity of the problem, 
we have {X\(t)Xi(t')) — (-Xi(O)-Xi(t = t' — t)), and the analytic expression of this autocorrelation function can be 
obtained via the Fourier transform of the spectrum of X\ Sx.iM = |x[^m + w]| 2 5j[u;] (where Sp[w] = ^MYyiksT 
is the thermal force spectral density): 

+00 



/+00 
-OO 



Similarly, noting that Jcj] = cu Sx,i[w]i one obtains the autocorrelation function of X^: 

- w^"-^-*" 1 )' (S7) 

where S is for the Dirac delta function. To determine the cross correlations (X1X1) and (X1X1), we write the quadra- 
ture Xi as the convolution between its impulse response x(t) = f_£> dwe l "'x[^M +uj] = 1 / (2MQ,M)Q{t)e (8 
denoting the Heaviside step function) and an effective thermal force Fh (with (Fth(i)Fth(t')) = Sp6(t — t')) [5]. 
After a tedious but simple calculation, we obtain: 

(S T ) 2 p-r M |t| 

(X 1 (t)X 1 (t'))(X l (t)X 1 (t')) = - -^ , 9 x— — . (S8) 



Combining Eqs. [S5j [S6j |7J and [S8j we finally find: 



l{e(tW)){fttW)) = (2 ^ f x^)e-»l'l. (S9) 

The autocorrelation of e remains to be calculated. To do so, we first note that (e(O)e(r)) = 2 ((X 2 (0)X 2 (t)) + (X 2 (0)) 2 ) 
The equipartition of the energy gives immediately the second term of the previous parenthesis, (X 2 (0)) 2 = 
(feaT/Mf^) 2 . The first one, that is the autocorrelation of X 2 , remains to be determined, which we do start- 
ing with expanding it as: 

(A 2 (0)X 2 (r)) = JJJJ dt 1 dt 2 dt 3 dt4x(0-ti)x(0-i 2 )x(r-t 3 )x(r-i 4 )(i ;i th(ii)i ;i th(i2) J F 1 th(i3)^th(i4)), 

[-oo,+oo] 4 

(S10) 

(Fth(ti), F t h(t 2 ), -Fth (£3), ^thC&O) being a multivariate Gaussian distribution, Wick's theorem applies, and the sta- 
tistical average of Eq |S 10| can be rewritten as: 

(Ft h (h)F th (t 2 )F th (t 3 )F th (U)) = (F th {tx)F th {t 2 ))(F th {t 3 )F ih (t 4 )) + (F th (ti)F th (t3)XF th (t 2 )F th (t 4 )) 

+<F th (ti)F th (t 4 ))<Ft h (t a )F th (t 3 )). (Sll) 

The previous equation shows that the fourth moment of the thermal force is a combination of products of its 
autocorrelation functions (F t h{ti)F t h(tj)) = 5p [f2 ~ £l m )5(ti — tj): 

(Fth(ti)Ft h (t 2 )Ft h (t 3 )Ft h (t4)) S&[fim] 2 {S(h - t 2 )6(t 3 - u) + 5{h - t 3 )S{t 2 - U) 

+5(t 1 -U)8{t 2 -t 3 )}. (S12) 



12 



Replacing Eq S 12 in Eq S 10 one obtains: 



T\2 



(Xf(0)Xf(r)) = (Sp 1 ) 



/+oo \ ■« / <.+oo 

dt lX 2 (*i) +2 / dt 
-OO / \J —oo 



X(h)x(r-ti) 



(S13) 



Using the expression of the impulse response x{t) given above, one finally obtains for the autocorrelation of e: 



(e(O)e(r)) 



( k B T 



M 



l + e 



-r M |t| 



(S14) 



Last, using that S<ptp[u] = & Eqs. S9 and S 14 lead to the very simple expression for the phase noise: 



(S15) 



The phase noise of a free-running, thermally driven mechanical resonator is hence entirely determined by and 
proportional to its dissipation rate, which is a well known result. Note however that importantly, this does not 
depend on the temperature. 



1.2 Feedback driven oscillator at thermal equilibrium 

It is straight to show that the expression |S 15| extends to the case of a feedback driven oscillator at thermal equi- 
librium: Indeed, the effect of linear dissipative feedback is to change the susceptibility to an effective susceptibility 
Xfb [u>] with the linewidth Tm being replaced by an effective linewidth Tfb which is an affine function of feedback 
gain. Noting that the thermal force is unchanged by feedback, and assuming that the feedback noise is negligible, 
the phase noise of the feedback driven resonator is obtained by deriving the above analysis with replacing Tm by 
Tfb, which leads to: 

CM = §r- ( S16 ) 

Thereby, positive dissipative strong feedback (r c ff <C Tm) enables considerable reduction of the phase noise of the 
resonator. While this result is somehow similar to the situation of an externally driven resonator as we will see 
in the next section, the phase noise decrease is limited to the stability of the external reference, while feedback is 
self-referenced. The experimental performance of the feedback technique will therefore mostly rely on the ability 
to design a low- noise feedback circuit. 



1.3 Externally driven resonator 



To do so, we start back with Eq. |S 2| 

We now turn to the case of an externally driven mechanical resonator. We further assume that the oscillator is 
resonantly driven using a coherent amplitude {X 2 ^) 1 / 2 > {k B T /M^) 1 / 2 . As a consequence, since Xi has been 
arbitrarily defined as being the phase reference (see Fig. |S1[ ), its (random) thermal component can be neglected 
Xth.i (Xi(t) = (^osc) 1//2 + -^th.l(i) — (^osc) 1 ^ 2 )- The motion phase imprecision resulting from the thermal 
fluctuations X t h,2 of the quadrature X2 is thereby given by: 

^"(A^' (S17) 

The phase noise associated to an externally driven state at thermal equilibrium therefore follows immediately from 
Eq. [STF] 

= yf§^, (S18) 

which is a well known result [TJ [5] : As increasing the external drive amplitud^] the thermal induced phase noise 
decreases. However, as noted above, the frequency stability limit then relies on that of the external reference. 
Reaching better stability performance will be thereby preferably possible using the positive feedback technique as 
announced above. 



1 We assume the mechanical motion to remain linear. 
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Figure S2: Phase and frequency stability of a mechanical resonator at thermodynamic equilibrium, (a) Schematic 
of the displacement spectrum of a thermally-driven oscillator (blue), feedback driven resonator with (X^) = (A t 2 h ), 
and externally driven resonator (brown) with (X% sc ) — (X^ h ). (b) Corresponding phase noise spectra, (c) Phase 
noise of an externally driven high-Q mechanical oscillator as a function of the Fourier frequency and mechanical 
damping-rate (normalized to an arbitrary damping unit Tq). (d) Phase noise of an externally driven high-Q 
mechanical oscillator as a function of the Fourier frequency and driving power. In both (c) and (d), the phase 
noise decreases towards darker regions, (c) and (d) show that increasing the external drive or mechanical quality 
factor lead to similar improvement for the phase noise performance, (e) Allan deviation of a thermally-driven 
oscillator (blue), feedback driven resonator with (X^) = (X^ h ), and externally driven resonator (brown) with 
(X% sc ) = (A t 2 h ). (f) Allan deviation for an externally driven resonator as a function of gate time and mechanical 
damping rate (normalized to an arbitrary damping unit Tq). The Allan deviation decreases towards the darker 
regions. Note the contours which tend to be parallel at higher damping rates, which reflects the fact that the 
frequency stability becomes independent of damping rate Tm at higher gate times. 



2 Thermodynamic limit of frequency stability 

In this section, we first establish the mathematical connection between phase noise and frequency stability. We 
then derive the analytical expression of the Allan deviation corresponding to the three cases treated above. 



2.1 Allan variance as a function of phase noise 

Here we establish the mathematical relationship between the Allan variance and the phase noise. For a given gate 
time r, the Allan variance <t(t) defines as: 




(S19) 
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where /m = ^m/^tt denotes the resonator's mechanical resonance frequency, expanding the squared parenthesis in 
the above equation, one obtains: 

a 2 (r) = -^Hmj^XM // dtdt' f(t)f(t') + Jj dtdt' f(t)f(t') - 2 JJ dtdt' f(t)f(t') 



(S20) 

where Dj lj2 = [(ji — 1)t, j'it] x [(j 2 — 1),J2t]. The latter equation can be rewritten using fixed integration's 
boundaries as follows: 

a 2 (r) = j — j-y E ( // dtd *'/(* + fcr )^' + kr)+ JJ dtd *f(t + (* - + (* - 1)^) 



fc=2 



'[0,r] 2 [0,r] ; 



2 yy didi'/(i + fcr)/(t' + (k - l)r) ) } . (S 21) 

[o.rp 

Swapping the statistical and temporal summations and using the ergodic assumption, one obtains: 

° 2 {r) = j^JJdtdt'({f(t)f(t'))-{f(t)f(t'-r))). (S22) 
M [0,tP 

To express the Allan variance in terms of frequency noise Sg [fi], we use the Wiener-Khintchin theorem which 
enables expressing the autocorrelation functions featured in the above expression as the inverse Fourier Transform 

+oo 

of the spectrum, (f(t)f(t')) = J dflSg [fl]e in ( t ~ t >. The latter equation reads therefore: 

— oo 

+oo 

<j 2 {t) = * r2 J J dtdt' J d£lS ff [0]e^'- f > (l - e- inr ). (S23) 

M [0,t] 2 -oo 

Fubini's theorem enables swapping the time and frequency integration to finally have: 

<7 2 (r) = _±^idQ^l(l-e-^)(l-cos(Or)). (S24) 



'M 

— oo 



Using that the imaginary part of the latter equation cancels, that the frequency spectrum is a pair function 
Sff[f2] = S'ff [ — S7] , and that the frequency noise S^ft] = Ss[ft]/fl 2 , the Allan variance finally reads: 



2 +°° 



«\r) = |(^_) /d^]si n *(^). (S25) 
o 

Though being similar to what can be found in Cleland et al [?, ?], this expression is a factor of it lower, which is of 
great importance: as an example, when the first one leads to the conclusion that a system is at the thermal limit, 
we find that it is actually still a factor > 3 above. 

2.2 Frequency stability of a mechanical resonator at thermodynamic equilibrium 

Using the result given by Eq. |S 25[ we are now able to determine the frequency stability associated to the cases 
treated in sections II. 11 11.21 and section 11.31 
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Free-running thermally driven oscillator Using Eqs. S 15 and S 25| one obtains for the frequency Allan 
deviation: 



(S26) 



with Qm = Hm/Fm denoting the mechanical quality factor of the resonator. Eq. S 26 shows that the stability 
increases with mechanical quality factor, gate time, and frequency: Therefore, besides being of high interest for 
purposes such as mass sensing |14[ |5]. high-frequency, high-Q mechanical resonators are extremely useful for absolute 
stability applications such as time and frequency control [TT] . 



Feedback driven oscillator at thermal equilibrium As suggested by the very similar mathematical expres- 
sion of the phase noise in the free-running and feedback driven configurations, the expression of the thermal limited 
frequency stability in the latter case takes the same form as Eq. 



S26 



a 2 (r) = 



1 



Qfb^MT ' 



(S27) 



with Qfb = f^M/rfb the effective mechanical quality factoij^] We note that Eq. S27 can be expressed differently, as 



a function of the feedback driven motion variance (A, 



th,fb/ 



4( 



r) 



(A 2 



th.fb/ 



1 



(S28) 



where we used that T^/Tm = (A t 2 h )/(A t 2 h fb ) 26 . This expression identifies to the one we derived in the next 
paragraph in the case of an externally driven oscillator at short gate times (r < 1/Tm), with the important difference 
that it applies at all gate times (in the limit r 3> 1/Om)- Again, the important difference with external driving is 
that this stability is no longer limited to the one of the external reference, which is found to be routinely on the 
order of a few ppb (10 -9 ) at second-scale for quartz oscillators. Reaching higher stability that become comparable 
to atomic clocks' (i.e. below 10~ 10 ) will therefore rather require the use of positive dissipative feedback. 



Externally driven resonator In the limit of a coherent drive 
driven resonator at thermal equilibrium is obtained using Eqs. 



S18 



> 
and 



S25 



frequency stability of an externally 
being given by [?]: 



with S'xx 



th r 

XX L 



F - 



(S29) 



We note that in MEMS-NEMS literature, this calculation is mostly achieved using a truncated expression of the 
phase noise, with the slowest noise component being neglected (S'xxM — ^xxl^ ^ Tm])- This is probably a 
technical legacy: MEMS technology first provided high-Q low-frequency oscillators, with very long mechanical 
coherence times 1/Tm, easily on the order of a few tens of seconds, along which technical limitations (e.g. thermal 
drifts, pressure fluctuations etc..) are the prominent sources of frequency noise. However, the recent development 
of high-frequency NEMS has led to low-noise, fast decay rate oscillators, the former truncation being henceforth 
obsolete. We thereby compute the integral in Eq. S 25 without further assumption, and find: 



knT 



Mni{xi, 



3 + cxp(-r M r) - 4 exp(-r M r/2) 



th/ 



2 \ 
osc/ 



x H(t). 



(S30) 



Eq |S 30| shows that the apparent frequency stability depends of the inverse ratio between the coherent drive and 
the thermal variance, along the thermal equivalent phase noise dependance. The function H{t) describes the gate- 
time dcpendancc, and features essentially 2 regimes. For short gate-times (t <C 1/Tm), we have H(t) ~ 1/Qm^MT, 
and the Allan variance identifies to the result obtained when neglecting the noise's slow components. In this case, 
one sees that a high Qm is favorable to lower measurement imprecision, along the usual optimization discussions. 
As already mentioned above, we also note that this asymptotic expression is identical to the effective frequency 
stability previously determined in the case of the feedback driven resonator, with (A 2 SC ) being replaced by (A^,), 
the use of feedback being though preferable as explained before. 



2 We assume the feedback mechanism to be purely dissipative, such that it does not induce any change in the mechanical resonance 
frequency. 
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On the other hand, for long gate times, Eq S 30 leads to H (t) ~ 3/i7|jT 2 : The thermally-limited Allan variance 
no longer depends on the mechanical quality factor, and decreases with the square power of r and with the 
fourth power of the mechanical resonance frequency. This 1/t 2 convergence behavior can be understood since 
averaging over gate-time greater than r ^> 1 /Tm means that two consecutive averaged realizations of the frequency 
measurement are independent due to the random nature of thermal noise, in contrast to the short gate-time case 
where these two measurements are related via the mechanical coherence time. This result also pleads for developing 
high frequency oscillators [51114) as highly accurate frequency references, even disregarding their mechanical quality 
factor, which should be kept high enough only for readout constraint^ 



3 The mechanical quality factor has to be such that the thermal noise can be resolved, since the frequency measurement imprecision 
will be due to detection background otherwise 
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